Grain boundary dynamics in stripe phases of non potential systems 
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^ . We describe numerical solutions of two non potential models of pattern formation in nonequi- 

librium systems to address the motion and decay of grain boundaries separating domains of stripe 



\ configurations of different orientations. We first address wavenumber selection because of the bound- 

ary, and possible decay modes when the periodicity of the stripe phases is different from the selected 
wavenumber for a stationary boundary. We discuss several decay modes including long wavelength 
undulations of the moving boundary as well as the formation of localized defects and their sub- 
sequent motion. We find three different regimes as a function of the distance to the stripe phase 
threshold and initial wavenumber, and then correlate these findings with domain morphology during 
Q ' domain coarsening in a large aspect ratio configuration. 
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I. INTRODUCTION 

"■^J , A stripe configuration is a particular example of a modulated configuration that exhibits structural periodicity along 
only one spatial direction. In thermodynamic equilibrium, these phases have been widely observed in a variety of 
systems [J, including, for example, microphase separated block copolymers @. Outside of equilibrium, stripe patterns 
i ^ i . are also widespread, including Rayleigh-Benard convection in fluids Q , or electroconvection in nematic liquid crystals 
When the lateral size of the system of interest is large compared with the wavelength of the modulation (large 
aspect ratio limit), a stripe configuration in either class of systems is usually characterized by a transient multidomain 
and defected configuration with only local order. Domain coarsening is believed to be controlled by the motion of 
existing topological defects in the configuration such as grain boundaries, dislocations, and disclinations. In this paper 
we focus on non potential (or no variational) model systems and investigate the processes of defect formation and 
\ motion near grain boundaries, and show that there exist several regimes in parameter space that are characterized by 
qualitatively different coarsening behavior. 
\Q ' A major topic of interest in domain coarsening of stripe patterns is the growth law for the characteristic linear scale 
of the configuration. This has been addressed both in theoretical studies [I, [1, 0, H, and in experiments in nematics 
d, [l(| and block copolymers [llj . This characteristic length scale R(t) is expected to obey a power law growth in 
time: R(t) ~ t x , with a value of the coarsening exponent x which ranges from 1/5 to 1/2 in the studies to date. The 
actual value appears to depend on e (the dimensionless distance from the onset of the stripe phase) [g, [ljj , thermal 
noise ||, and the type of length scale under study @, Also, domain morphologies as well as the value of x have 
been found to be different in potential systems (in which the evolution is solely driven by minimization of a potential) 
than in the non potential case [fjj . Numerical studies of the Swift- Hohenberg model of Rayleigh-Benard convection (a 
. potential system) suggest that x = 1/3 in the limit e < 1 as a result of grain boundary motion]?]]. As e increases, the 
value of x is seen to decrease from 1/3, fact that was argued to be a manifestation of defect pinning [12] . On the other 
hand, the generalized Swift-Hohenberg model has been studied as an example of a non potential model @, and it has 
been shown that the evolution of stripe patterns is asymptotically dominated by the motion of isolated dislocations. 
Furthermore, recent results from a direct solution of the Boussinesq equations in a Rayleigh-Benard configuration 
suggest that multiple length scales should be introduced in the case of non potential systems @ . 

Despite a significant body of research on domain coarsening of modulated patterns, the mechanisms responsible for 
the observed properties and the values of the coarsening exponents, especially for non potential systems, are not yet 
understood. We report in this paper our analysis of the relation between the selected wavenumber at grain boundaries 
and the critical wavenumber, and consequences on both early time defect formation and asymptotic coarsening. Our 
starting point is a planar 90° grain boundary separating two stripe domains of mutually perpendicular orientations, 
as shown in Fig. [T] We determine the selected wavevector for a stationary boundary configuration q s in two different 
non potential models, and show that for other wavenumbers the mode of decay of the configuration can be classified 
according to its wavenumber and the dimensionless distance from threshold e. For small e, an initially planar grain 
boundary separating two domains, one of which at least has q ^ q s , is seen to propagate, either as a planar front, 
or by developing long wavelength undulations, depending on the value of q. For larger values of e we find that the 



FIG. 1: A 90° tilt grain boundary (GB) configuration of a stripe phase, with wave numbers q x and q y for two domains A and 
B of mutually perpendicular stripe orientations. 



grain boundary configuration decays through the formation of disclinations around the boundary and their subsequent 
motion. We then correlate our findings with the observed domain coarsening morphology of configurations obtained 
from initially uniform but unstable states. The rationale for the comparison is the observation that in an extended 
system an early time transient configuration following the decay of an unstable state will comprise a distribution 
of locally ordered domains with wavenumber approximately equal to the critical wavenumber for instability go, but 
with different orientations. Depending on whether go is close to q s or not, the resulting boundaries will be locally 
close to being stationary, or quickly decay leaving behind different types of defects. The ensuing domain coarsening 
is shown to be different depending on the (near) stability of grain boundaries that spontaneously form in the early 
stages following the decay of the initially unstable configuration. 

Our study is based on a numerical solution of two different non potential extensions of the Swift- Hohcnberg model 
of Rayleigh-Benard convection as described in Sec. [Til One of them is made non potential on account of the 
specific choice of nonlinearity, whereas the other allows for mean flows representing the effect of vertical vorticity 
encountered in Rayleigh-Benard convection in small Prandtl number fluids [l3j . They are described in Sections 11111 
and IIVI respectively. Although wave number selection for a stationary grain boundary is different in both models, 
similar properties concerning grain boundary motion and decay have been found. They correlate with two types of 
domain coarsening behavior, either grain boundary dominated or disclination/dislocation dominated, as deduced from 
the results on grain boundary decay that will be discussed in Sec. [V] 



II. NON POTENTIAL MODEL EQUATIONS 



The model systems we have investigated are two dimensional gene ralized Swift-Hohenberg equations which have 
been used as a phenomenological description of stripe patterns (S Il3|: 

dfi/i + U- Vt/j = exjj - (V 2 + qlfi) - gifi 3 + 3(1 - .g)|WfV 2 i/;, (1) 

where ip{x, y, t) is a dimensionless real order parameter field, and the control parameter e is the distance from the 
threshold for instability of the solution tp — (at e = 0). The wavenumber go (go = 1 in the dimensionless units that 
we are using) is the critical wave number for instability. An advection term with a drift velocity U(x, y, t) = (U x , U y ) 
has been added to Eq. ([TJ to introduce mean flow effects resulting from the coupling to vertical vorticity in the fluid 
[T3L [lil]. The velocity U(x, y, t) is given in terms of the vorticity potential £(x, y, t): 

U = V x (Cz) = (d y C, -d m Q , (2) 

which satisfies, 

[8t - <j(V 2 - c 2 )] V 2 C = g m [V(V 2 ^ x VV] ■ z, (3) 

with <7 being the Prandtl number of the fluid, c 2 introduced to phenomenologically model the effect of no slip boundary 
conditions at the top and bottom plates of the convection cell |14l|. and the coupling parameter g m chosen to be a 
decreasing function of the Prandtl number. Following Ref. [lj a Gaussian filtering operator F 7 is introduced on 
the r.h.s. of Eq. ([3]) to reduce high wavevector contributions to the vertical vorticity. In Fourier space, acts by 
multiplying the Fourier transform of the filtered function by a factor 
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where q 2 = q 2 + q 2 , with wavevector components q x and q y , and 7 represents the filtering radius. 

Different types of model equations for describing stripe patterns have been included in Eqs. ((T])-©, depending on 
the values of g and g m . For g = 1 and g m = (so that U = 0), the nonlinearity in the model is given by the cubic term 
tp 3 only, leading to the original Swift-Hohcnberg equation [15| , which is a potential model. A non potential model 
without mean flows follows from g = g m = (referred to as model 1 in this paper), which includes a nonlinear term of 
the form | Vip\ 2 V 2 ip. For g — 1 and g m ^ (and hence U ^ 0), we obtain the Swift-Hohenberg model supplemented 
by mean flows, also a non potential system (model 2). 

The initial configuration considered in the first part of our study is a 90° tilt boundary separating two stripe 
domains, as shown in Fig. [TJ A stationary state is known to exist [l6| if unique values of the wavenumbers q x (for the 
stripes parallel to the grain boundary in domain A) and q y (for the stripes perpendicular to the boundary in domain 
B) are selected in both domains. For the original Swift-Hohcnberg model, these wavenumbers are q x = q y = qo. In 
our numerical study, we consider instead a symmetric pair of grain boundaries along the x direction so that periodic 
boundary conditions can be adopted. The pseudospectral algorithm described in Ref. [13 has been used to numerically 
solve the model equations, with a time step At = 0.2 and a grid spacing Ax = Ay = Ao/8 (i.e., 8 grid points per 
stripe wavelength Ao = 2ir/qo). Most of our results shown below correspond to a system size 512 x 512 grid nodes, 
with spot checks with larger sizes (2048 x 2048) to verify our results. 



III. WAVEVECTOR SELECTION AND GRAIN BOUNDARY DECAY - MODEL 1 



A. Wavenumber selection 



Model 1 includes the nonlinear term |V?/>| 2 V 2 ?/> but does not have mean flow. The selected wave numbers (q x and 
q^) for a stationary grain boundary configuration have been already calculated by Tesauro and Cross (l6| . Their 
analysis of the appropriate amplitude equations for this case yields 

€ = *> - «3 ( 4 ) 

up to order e for the selected wavenumber in domain A, and q y = qo up to order e 3 / 4 for domain B. They also 
conducted a limited numerical study of this model to suggest that q y ~ go — 0.1 2e (although this result is based on 
calculations involving only two different values of e). 

We have numerically calculated the selected wave numbers for the same configuration but with much larger system 
sizes and longer integration times. Our results for g* are in agreement with Eq. (J4|, but the calculated values of 
q^ are different from those of previous studies, as presented in Figs. [5] and [31 The deviation of q y ' from the critical 
wavenumber go is seen to obey two different power laws in e. For small e < £qi with eo around 0.2, we find (Fig. [3J 

q s y ~ go - 0.23e, (5) 

while for larger e we obtain, 

q s y ~ q - 0.525e 3 / 2 . (6) 

Note that although both results differ from the numerical calculations of Ref. [TH, they are not inconsistent with an 
analysis based on amplitude equations. 



B. Decay of the grain boundary configuration 

Given the results of Eqs. (2])-©, and also our numerical results in Fig. [51 it is clear that the wave number selected 
by a grain boundary is different from go, the most unstable mode near onset. Thus the locally dominant wavenumber 
during the early stages of coarsening could be significantly different from (g*,g*), and therefore it appears to be of 
value to study the decay of grain boundaries under those conditions. Of course, grain boundaries would not be fully 
formed at those early stages, but we hope that studying unstable trajectories away from well defined grain boundary 
configurations will illustrate differences in observed morphologies between potential and non potential models. 

We start by noting that in the configuration which we consider there is no phase conservation in the x direction. 
Hence if q x ^ q a x is chosen initially, the bulk wave number of domain A can readjust to g* through lamella formation 
or destruction near the boundary [r|. We will not consider this case. More complicated phenomena can occur for 
% ^ Qy' an d this is the focus of our study: In the calculations that follow we always keep q x = q x in domain A, while 
we explore a range of values of q y in domain B. 
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FIG. 2: Stability diagram for the non potential model 1 as a function of e and q y (the wave number of the stripes perpendicular 
to the grain boundary). The plus symbols correspond to the unique wavevector (q x ,q y ) = (gj = qo — e/(2qo), Qy) for a 
stationary grain boundary configuration. Also shown (as dot-dashed lines) are stability boundaries for bulk Eckhaus and 
zig-zag instabilities. Three regimes for grain boundary decay are indicated, and also the critical wave number qo for linear 
instability. 




FIG. 3: Wave number q y for a stationary grain boundary (model 1) as a function of e. The figure shows two different types of 
dependency on e: go — q y cx e (for small e) and oc e 3 / 2 (for larger e). 




FIG. 4: Grain boundary decay in the three regimes shown in Fig. [2] for model 1, with q x = qj. The direction of boundary 
motion is indicated by arrows, (a) Regime I: q y < q y for all e; the example shown here corresponds to q y = 0.85, e = 0.3 and 
at t — 1200. (b) Regime II: q y > q y and e < e-r; we have chosen q y = qo = 1, e = 0.04, and t = 10 4 . (c) Regime III: q y > q y 
and e > ex; the configuration shown here corresponds to q y = qo = 1, e = 0.5, and t = 100. 



(a) t=0 



(b) t=500 



FIG. 5: Grain boundary decay (for model 1) due to a slow transverse modulation of domain A in an initial stationary 
configuration (with q x = q x and q y = q y ), with e = 0.5, the modulation wavenumber Q — q x /16, and the magnitude of the 
phase modulation 5 — Xq. 

As indicated in Fig. [21 we find that our results can be qualitatively classified according to the values of q y and e. 
In regime I defined by q y < q y , the two grain boundaries move from domain A towards domain B until they merge. 
Both boundaries remain approximately planar during the process. Such a motion is found for all the values of e 
investigated, with a particular example of e = 0.3 and q y = 0.85 shown in Fig. 2^. 

On the other hand, a longitudinal distortion of the stripes in domain B around the grain boundary region is 
observed if q y > q y . In addition, the evolution following the appearance of the initial distortion is qualitatively 
different depending on the value of e. In regime II of Fig. [5] with e < ey (where the value of ey is around 0.1 to 
0.2 as estimated from our numerical solutions), long range undulations of diffuse grain boundaries are observed, as 
shown in Fig. [4b. The stripes of domain B are curved around the boundaries, the latter moving towards each other 
at the expense of domain B until they finally merge. Despite the undulations, the grain boundary configuration in 
regime II is preserved during the process. In regime III on the other hand, the boundary is sharp, and the longitudinal 
distortion leads to the formation of convex disclinations in the boundary region (Fig. for e > ey). The boundary 
itself also moves towards the bulk of domain B as indicated by the arrows in the figure, resulting again in the eventual 
disappearance of the grain boundary configuration. 

The longitudinal distortions that we observe in both regimes, including the formation of an array of disclinations 
in regime III, are due to wavenumber mismatch between the initial value q y and the stationary value q y < q y . In 
regime II (Figs. Hb), a periodic longitudinal modulation of stripes B appears, reminiscent of an Eckhaus distortion 
that promotes wavenumber reduction. Since the velocity of the boundary depends on the local wavenumber of stripes 
B ahead of it, a modulation of the boundary position also results from this wavenumber variation. In regime III, on 
the other hand, the boundary is sharp, and cross rolls, dislocations and disclinations form at the edges of domain B 
leading also to wavenumber reduction. The bulk of domains A and B, however, remain largely unperturbed. 

In order to verify that this spontaneous distortion originates in the boundary region, we have studied the decay of 
the configuration shown in Fig. [5^ in which slow transverse modulations of wavenumber Q <C q x and phase magnitude 
8 have been added to the initial condition of q x — q^. in domain A. In this case, q y = q y in the bulk of domain B but not 
at the grain boundaries. For large e (regime III), the grain boundaries decay through the formation of disclinations, 
as shown in Fig. [Sb. The same decay was shown in Fig. 3J; corresponding to the case q y ^ q^ in the whole domain B. 

A similar analysis can be carried out in regime II (small e), with the same conclusion. We note, however, an 
interesting aspect of the decay as compared to potential models. The grain boundary moves from the modulated 
domain A towards the undistorted domain B, exactly the opposite of what is expected for a potential model given 
our choice of wavenumbers. With the consideration that wavenumbers q x and q y in the bulk of both domains are 
equal to the stationary values and q y of the corresponding potential model, the extra energy that is stored in the 
modulated stripes in domain A would lead to boundary motion in the opposite direction 0]. 

IV. WAVEVECTOR SELECTION AND GRAIN BOUNDARY DECAY - MODEL 2 

The second non potential model studied, as described in Sec. HH includes the so called mean flows. They are known 
to have a significant effect on the stability of a stripe configuration [l3[ , on defect dynamics [H, Ell > an d allow the 
appearance of spiral defect chaos [2(| • We consider here the decay of the grain boundary configuration studied above. 
We have chosen model parameter values that are typical for fluid systems in Eqs. ([I])-©: <r = 1, c 2 = 2, g m = 10, 
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8 < e x e > e T 




(a) Regime II (b) Regime III 

FIG. 6: Grain boundary decay for model 2 as a result of mean flows, with g m — 10, a — 1, <? — 2, and q x = q% — qo. (a) 
Regime II: q y > q y and e < €t', the parameters used in this case are q y = 1.02, e — 0.04, and t — 4000. (b) Regime III: q y > q y 
and e > er; the configuration shown corresponds to q y = 1.02, e = 0.3, and t = 600. 

and 7 = Ao/2 for filtering. 

We have first verified numerically that the unique wavenumbers for a stationary grain boundary are = qt = qo, 
the same as in the potential Swift-Hohenberg model. This result follows from the observation that the lowest order 
contribution of mean flow to the amplitude equations of a grain boundary is C(e 5 / 4 ) [lj]. 

The known stability diagram of a stripe configuration is similar to that shown in Fig. [5] for model 1, except that 
Qy = Qo = 1, an d the type and detailed location of the stability boundaries are different (e.g., the Eckhaus boundary 
in Fig. [5] should be replaced by the skewed- varicose boundary [138). I n terms of our study of the decay of a grain 
boundary, we also find in this model three different regimes which are qualitatively similar to those of model 1. Rigid 
motion of the grain boundary is observed when q y < for all e, while a longitudinal distortion of the stripes in 
domain B appears when q y > q y . We can also qualitatively distinguish between regimes II (e < ex) and III (e > er), 
with two typical examples given in Fig. El 

As was the case in model 1, we observe long wavelength distortions of the stripes in domain B in regime II (Fig. 
[S^i), whereas regime III leads to the formation and subsequent motion of disclinations, and with them, the decay of 
the grain boundary configuration (Fig. [6)d) . Disclinations also appear as a result of local deviation of wavenumbers 
from q^ in the grain boundary region. This has been verified by adding slow transverse modulations on A stripes 
while keeping q = q s in the bulk of both domains, with results that are qualitatively similar to those of Fig. 

The driving force behind the appearance of longitudinal distortions, and with them either long wavelength undula- 
tions or the formation of disclinations is again local wavenumber reduction around the boundary. This is analogous 
to what is observed in model 1, and has been already described in detail in Sec. IIII Bl However, once the undulations 
have been formed or the defects nucleated, mean flows are effective in their subsequent motion. This can be seen from 
the spatial distribution of the vorticity potential £, or the corresponding drift velocity U obtained in our numerical 
solutions. In particular, we note that mean flows play a key role in the motion of disclinations in regime III, a regime 
in which disclinations become pinned in their absence {g m — > 0). 

V. DISCUSSION AND CONCLUSIONS 

The results described on the decay of a grain boundary configuration in non potential models shows the existence 
of a range of parameters in which the decay leads to the formation of defects (mostly dislocations and disclinations) 
in the boundary region. This occurs for a configuration in which the stripes parallel to the boundary plane have 
the selected wavenumber q x = q s x , whereas the stripes perpendicular to the boundary are under slight compression 
q y > q y . Also, defects are only seen above a certain threshold e > ex- Below the threshold, undulations and motion 
of the boundary follow, but no additional defects are produced. 

These results can help interpret qualitative features of domain coarsening behavior, both in potential and non 
potential models. In the case of non potential model 1 (with a | V"0| 2 V 2, nonlinearity) , the stationary wave numbers 
q x and q y are different from qo, the most unstable wave number (see Eqs. (JU)-© and Fig. Following the decay 
of an initially disordered configuration (tp = 0), and the appearance of locally ordered domains with wavenumber 
approximately equal to qo, grain boundaries would be expected to decay according to the mechanism described for 
regimes II and III depending on the value of e. Whether coarsening takes place in the regime in which the decay 
of boundaries leads to the formation of dislocations and disclinations (regime III), or to their motion (regime II) 
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(a) 8 < 8 X (b) 8 > 8 X 

FIG. 7: Multidomain configurations for model 1 showing different qualitative coarsening behavior for (a) e < er (with e = 0.04 
and at t = 1000) and (b) e > e T (with e = 0.5 and at £ = 500). 




(a) 8 < 8 X (b) 8 > 8 X 

FIG. 8: Multidomain configurations of stripe phase in non potential model 2 (with mean flows), for (a) e < et (e = 0.04 and 
£ = 1000) and (b) e > er (e = 0.5 and £ = 1000). Other parameters are the same as those in Fig. g m = 10, a = 1, and 
c 2 = 2. 

one might anticipate qualitatively different late time domain coarsening. Figure [7] shows the results of a numerical 
integration of model 1 from a initial random configuration of xjj for two particular values of e, both below and above 
ex- For e < £t (Fig. EK)> grain boundaries are visible in the late time configuration as expected from the results shown 
in Figs. [4Jd. Qualitatively different behavior is observed for e > ex due to the early stage formation of dislocations or 
disclinations around the emerging grain boundaries. Our results in this case are shown in Fig. [TJd, with the resulting 
disclination/dislocation dominated morphology. These results are consistent with earlier observations in Ref. [6|, and 
more recent research on the Boussinesq equations in Ref. [^. The average wave number of the evolving configuration 
changes from go at early times following the instability, towards a long time value that is closer to the stationary wave 
number of isolated dislocations. One anticipates that domain coarsening will be qualitatively different from the case 
of small e, as suggested by the results of Ref. [a showing a complex scaling behavior at large e. 

In the case of model 2 which includes the coupling to mean flows, the selected wavenumbers q% and q y approximately 
coincide with the critical wavenumber for instability qo, and therefore the effects described above will be weaker, but 
not completely absent as local deviations between the stripe wavenumber and the selected wavenumber still occur in 
regions in which the stripes are curved. Near threshold (e < €t), the morphology of a transient stripe configuration 
is shown in Fig. [5^, obtained by numerical solution of the governing equations from a random initial condition. A 
stripe pattern qualitatively similar to both model 1 and the potential Swift-Hohenberg model is found, characterized 
by differently oriented domains separated by 90° grain boundaries. On the other hand, when e exceeds a certain 
threshold er, grain boundaries are expected to be absent as described in Sec. IIV1 The subsequent evolution of the 
system in this case of is much more complicated resulting in persistent dynamics (spiral defect chaos (201). as shown 
in Fig. 03. 
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In both non potential cases discussed in this paper we would expect qualitatively different domain coarsening 
behavior depending on the interplay between the critical wavenumber for instability and the selected wavenumber at 
a grain boundary. This is in contrast with the potential Swift- Hohenberg model in which grain boundary dynamics 
is known to dominate for small e until pinning effects at finite e effectively arrest motion y§, [H[ • In this model, the 
selected wavenumbers for a stationary grain boundary configuration are q% = q y = Qo [l6j . with q = qo also marking 
the stability boundary of the bulk zig-zag instability. Therefore a diagram analogous to that of Fig. [2] would only 
have regimes II and III characterizing grain boundary motion. In both q y > q^, longitudinal distortions around the 
boundary can also be identified as a result of local wavenumber reduction towards q y = qo . In analogy with the non 
potential cases, in regime II (with e < ex) diffuse boundaries with long range undulations of B stripes are observed. 
For e > €t, the formation of convex disclinations around the boundary can also be seen in numerical calculations. 
However, despite their emergence these disclinations become pinned along with the grain boundaries, mainly due to 
the nonadiabatic effects arising from the coupling between fast scales of base stripe patterns and slow scales of the 
associated envelopes (l2| . This results in a glassy configuration with very slow dynamics. Interestingly, this pinning 
phenomenon is only observed in the potential Swift-Hohcnberg model, but not in the two non potential models studied 
here. In this latter case, grain boundaries quickly decay at large e leaving behind a domain morphology dominated 
by disclinations and dislocations (Fig. [TJd), or spiral defects (Fig. [5}d). 

In summary, we have studied wavenumber selection, motion, and decay of grain boundaries in a stripe configuration 
based on a pair of two-dimensional generalized Swift-Hohenberg model equations that are non potential in nature. 
Three different regimes of grain boundary dynamics have been identified, with longitudinal distortions emerging 
around the grain boundary and manifesting themselves either as long range undulations of the moving boundary or in 
the formation and subsequent motion of defects such as convex disclinations. When an extended system is taken from 
an initial state below threshold to e > 0, the initial decay of the unstable state will lead to the formation of locally 
ordered domains with a distribution of orientations. Domain interfaces will form in accordance with the analysis given 
above, and therefore qualitative differences in domain coarsening are expected depending on the value of e. 
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